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We study the dynamics of a tracer particle (TP) on a comb lattice populated by randomly moving 
hard-core particles in the dense limit. We first consider the case where the TP is constrained to move 
on the backbone of the comb only, and, in the limit of high density of particles, we present exact 
analytical results for the cumulants of the TP position, showing a subdiffusive behavior ~ At 
longer times, a second regime is observed, where standard diffusion is recovered, with a surprising 
non analytical dependence of the diffusion coefficient on the particle density. When the TP is 
allowed to visit the teeth of the comb, based on a mean-field-like Continuous Time Random Walk 
description, we unveil a rich and complex scenario, with several successive subdiffusive regimes, 
resulting from the coupling between the inhomogeneous comb geometry and particle interactions. 
Remarkably, the presence of hard-core interactions speeds up the TP motion along the backbone of 
the structure in all regimes. 


Subdiffusive motion of tracer particles in crowded me¬ 
dia, e.g. biological cells, is widespread. Among the possi¬ 
ble microscopic scenarios leading to this sublinear growth 
with time of the mean square displacement (MSD), the 
existence of geometric constraints related to the complex¬ 
ity of the environment plays an important role [1, 2]. In 
this context, the comb model (see Fig. 1), in which a sin¬ 
gle particle moves on a two-dimensional space with the 
constraint that steps in the x direction are only allowed 
when the y coordinate of the particle positions is zero, 
has attracted considerable attention because of its sim¬ 
plicity and ability to reproduce subdiffusive behaviors of 
disordered systems [3]. 

Comb-like structures have indeed been introduced as 
a first step to model diffusion in more complicated frac¬ 
tal structures like percolation clusters, the backbone and 
teeth of the comb representing the quasilinear structure 
and dangling ends of percolation clusters [4]. The par¬ 
ticle can spend a long time exploring a tooth, which re¬ 
sults in a subdiffusive motion along the backbone with 
{x‘^{t)) oc with a = 1/2. Since, numerous results 
have been obtained for this model [5-12], including the 
determination of the occupation time statistics [13], of 
mean first-passage times between two nodes of a finite 
comb [14] or the case of fractional Brownian walks on 
comb-like structures [15]. 

In parallel, the comb model has been invoked to 
account for transport in real systems like spiny den¬ 
drites [10], diffusion of cold atoms [16] and mainly dif¬ 
fusion in crowded media like cells [17]. However, all ex¬ 
isting studies have focused on single-particle diffusion, 
and interactions between particles have up to now been 
completely left aside. As an elementary model for dif¬ 
fusion of particles under short-range repulsive forces, we 
consider here excluded-volume interactions (EVIs) and 
focus on their impact on tracer dynamics on comb-like 
structures. 

From a theoretical point of view, lattice systems of 
interacting particles represent a prototypical model in 
statistical physics that has generated a huge number of 
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FIG. 1. Geometry of the system. The a;-axis is designated 
as the backbone of the lattice, whereas the orthogonal lines 
are designated as the teeth. Jump rules of the particles in the 
case when the tracer particle (in red) is constrained to move 
on the backbone are given. 


works both in the physical [18, 19] and mathematical lit¬ 
erature [20]. The effect of EVIs in homogeneous lattices 
is well known [3]. In dimension d > 2, tracer diffusion 
has been shown to remain normal, with a non trivial 
diffusion coefficient resulting from many-body interac¬ 
tions and well approximated by the Nakazato-Kitahara 
approach [21]. In a “single-file” geometry, where particles 
cannot bypass each other, the impact of EVIs is stronger 
and results in a subdiffusive behavior {x^{t)) oc with 
/3 = 1/2 [22-27]. In this context, determining the effect of 
EVIs on systems with geometrical constraints appears an 
important question which does not seem to have received 
much attention. Notable exceptions are [28], where two 
particles only are involved, [29], which involves single-file 
motion with waiting times and thus does not consider 
explicitly geometric constraints, and [30], where tracer 
diffusion on DLA clusters was studied numerically and 
found to be not modified by the presence of EVIs. In 
this Letter, we show that, in contrast, EVIs deeply mod¬ 
ify tracer diffusion on comb-like structures. Focusing on 













2 


the high density limit, we show analytically that, due to 
a subtle interplay between the inhomogeneous comb ge¬ 
ometry and crowding effects, the dynamics displays sev¬ 
eral regimes of anomalous diffusion. We find in particu¬ 
lar that, surprisingly, the presence of EVIs can speed up 
tracer diffusion along the backbone of the structure. 

Model- We consider the 2-dimensional comb C 2 , 
which is a subgraph of 1? obtained by removing all the 
edges parallel to the a;-axis, except from the x-axis it¬ 
self. This lattice is populated by N hard-core particles 
with average density p = N/V, where V is the number 
of sites. Each particle performs a symmetric nearest- 
neighbor random walk, constrained by hard-core inter¬ 
actions. We add a tracer particle (TP) at the origin, 
performing a symmetric nearest-neighbor random walk, 
and focus on its dynamics in the dense limit, where the 
vacancy density po = 1 — p 1. In this limit, it is more 
convenient to describe the vacancy dynamics instead of 
describing the dynamics of all the particles. We assume 
here that, at each time step, each vacancy exchanges its 
position with one of the neighboring particles, with jump 
probabilities that depend on the position on the lattice; 
see the Supplemental Material (SM) [31] for the explicit 
definition of those evolution rules. 

Case of a TP restricted to the backbone- We first as¬ 
sume that the TP (and only the TP) is constrained to 
move on the backbone. This particular case is impor¬ 
tant for several reasons, (i) It mimics the case where the 
tracer is different from the bath particles, and is not able 
to visit the teeth of the comb, (ii) At the theoretical 
level, it appears as an extension of the famous single-file 
geometry in which, due do the possibility for the bath 
particles to visit the teeth, the particles can bypass each 
other. An interesting question is to know if the dynamics 
is still anomalous in this case, and if so, with which expo¬ 
nent; (iii) finally, as shown below, solving this auxiliary 
problem will allow us to determine the dynamics of the 
TP in the general case where the TP can access the teeth 
of the comb. 

Let Xt be the random variable characterizing the po¬ 
sition of the TP along the backbone at time t. We aim 
at computing the cumulants of arbitrary order n of this 
variable, denoted by K^”^(t). These quantities are gener¬ 
ated by the cumulant generating function (CGF) 'l't(fc) = 

ln(e®*^*) = In Pt{k) , where Pt{k) = Pt{X) is 

the Fourier transform of the probability Pt{X) to find the 
TP at position X at time t. Following the method devel¬ 
oped in [32] and recently used to study driven diffusion 
in one-dimensional geometries [33], we first consider the 
case where there is a single vacancy on the lattice. Let 
P^^\x\Z) be the probability to find the TP at position 
X at time t knowing that the vacancy started from site 
Z. Summing over all the passages of the vacancy to the 


TP location, one gets: 


P^^\x\Z)=5xp I 1-Y.FMZ) 
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where Ft{Q\Z) is the probability for the vacancy to reach 
the origin for the first time at time t, knowing that it 
started from site Z, and ei stands for the unit vector in 
the X direction. The first term in the right-hand side of 
Eq. (I) represents the event that at time t, the TP has 
not been visited by any vacancy, while the second one re¬ 
sults from a partition both on the number p of visits and 
waiting times mi between visits of the TP by the vacancy. 
Computing the generating function associated with this 
propagator p±i(Ar; 5) = P^^^X] ± ei; ^), where de¬ 
notes the discrete Laplace transform 
and noticing that for symmetry reasons F(0|ei;^) = 
^(0| -ei;^) = Fi, one gets 


P±iiX-A) = 


'Jx.o)! — -Fi) -I- (5j(:,±iFi(I — Fi) 


We then study the case where the concentration of va¬ 
cancies on the lattice po is finite but very small. For clar¬ 
ity, we first assume that the lattice has a finite number of 
sites N, and that it is populated by M vacancies, so that 
M = poN. Consequently, the probability Pt{X\{Zj}) to 
find the TP at position A as a result of its interactions 
with M vacancies initially located at sites Zi, • • • , Zm, 
is given by 


Ptix\{z,}) = sx,Y,+-+Y^Ptm}\{Zj}), 

Z 1 ,Zm 

( 3 ) 

where Pt{{Yj}\{Zj}) is the conditional probability that 
during the time interval t the TP moved of a distance Yj 
due to its interactions with the j-th vacancy. To lead¬ 
ing order in po, the vacancies contribute independently 
to the displacement of the TP, so that in Fourier vari¬ 
able, pj^'^ (fc) = (fc) j where (fc) is the Fourier 

transform of the probability distribution to find the TP 
at position X at time t, knowing that there are j vacan¬ 
cies on the lattice, and averaged over the initial position 
of the vacancies, which is assumed to be uniform. As 
shown below, the choice of the initial distribution of the 
vacancies may actually have a dramatic effect on the be¬ 
havior of the TP. 

Finally, in the thermodynamic limit in which A, M —)■ 
00 with fixed po = M/N, and using Eqs. (2) we get the 
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FIG. 2. (Color online) Variance rescaled by po computed 
from Monte-Carlo simulations of vacancy dynamics. The full 
line represents the analytical result (5). Inset: case where all 
vacancies are initially placed on the backbone (the full line 
represents Eq. (12.)) 


Fourier Laplace transform of the CGF 

~ - fu? 9, (l-cosfc), (4) 

Po^O + 

where we defined H{^) = E^-oo ^ 2/i 0: 

and where we used the symmetry relation 
Ez^o^*(0|ei|Z) = Ez^o^*(0| - ei|Z). Conse¬ 
quently, the determination of the CGF amounts to the 
calculation of the quantity H{^) and the first passage 
density which are given in the SM [31]. 

Expanding lE'(fc; ^) in powers of k from Eq. (4), focusing 
on the large time limit ^ 1“ and using a Tauberian 

theorem [34] we then get the following exact expression 
for the cumulants (see SM [31] for details) 

limdffW 1 (3/4 ,5) 

Several comments are in order, (i) Eirst, one notices that 
all cumulants of the same parity are equal. This indicates 
that the probability distribution Pt{X) is a Skellam dis¬ 
tribution [35], originally defined as the p.d.f. associated 
with the difference of two Poissonian random variables. 
In particular, the rescaled position is asymptot¬ 

ically normally distributed, (ii) The exponent 3/4 is in¬ 
termediate between a “normal” diffusion exponent, and 
the single-file diffusion exponent 1/2. (iii) While the 
analytical predictions of the cumulants are successfully 
compared to results obtained from Monte Carlo simula¬ 
tions at intermediate times (see Fig. 2), at long times 
a crossover towards a standard diffusive behavior is ob¬ 
served. 

We present below a theoretical argument that accounts 
for this intriguing behavior. The key point underlying 
this crossover is that the above analytical results are de¬ 
rived by taking the limit po 0 before the long time 


limit. In this case, to leading order, the TP does not move 
before being reached by a given vacancy. In fact, the TP 
actually diffuses due to its interactions with the other va¬ 
cancies. Because of this effective diffusion of the TP, each 
vacancy experiences an additional symmetric jump prob¬ 
ability in the x direction, denoted by D{pq), even when 
they are on the teeth. Thus, in the reference frame of the 
TP, the vacancies are allowed to jump from one tooth to 
another (with a probability that vanishes as po —> 0) and 
their motion is now effectively two-dimensional. Quali¬ 
tatively, we are facing a two-dimensional situation, and 
regular diffusion is thus expected at large times (see intro¬ 
duction). Quantitatively, the approach developed previ¬ 
ously can be extended to calculate the variance of Xt, 
in this case where the vacancies perform 2D random 
walks, but can reach the TP only from the backbone 
(see SM [31]), yielding 




(?) = -2po- 


S(e,Po) 


- - I 


F*-l 


. F* 


F* + 1- Ft 




( 6 ) 


where we defined F^i = F*(Ojeije±i;^,po) and 

F‘{^,Po) = Ez5^o^*(0|ei|^;^,Po), with F* {0\ei\Z-, po) 
the probability for a vacancy to reach the origin for the 
first time at time t knowing that it was at site ei at time 
t — 1 and that it started from site Z. Relying on renewal- 
type equations, the conditional first-passage time densi¬ 
ties F* are related to the propagators of the vacancies 
random walk, which are themselves evaluated with an 
extension of the method presented in [36] to treat diffu¬ 
sion on inhomogeneous lattices (see SM [31]). It is finally 
found that 


■^-1 - — VD{po) In 


Ty^D{po) 


D{po) In Ej 


Ff ~ 1 - 


WF>{po 


+ ... (7) 


( 8 ) 


where the symbol ~ stands for the long-time limit (^ —> 
1) for a fixed value of po- As expected, Fti —t 0. Next, 

po-70 

E(^,po) is conveniently estimated in SM [31] by relying 
on the continuous-space description given in [5]: 


m,po) 




( 9 ) 


Plugging Eqs. (7), (8) and (9) into (6), and using a 
Tauberian theorem, we finally obtain 


lim 

t—^OO 


t 


oc po\/D(po)ln 

Po-i-O 


1 

mp^)' 


( 10 ) 


This equation defines the effective diffusion coefficient 
D{po) = liint_>.oo(K^^^(t)/2t) self-consistently when po —>■ 
0, and finally yields the following expression of the vari¬ 
ance in the ultimate regime: 


lim 
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( 11 ) 
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These results thus show that the limits po —>■ 0 and 
t —>■ (X) do not commute leading to an ultimate diffu¬ 
sive behavior [37]. However, due to a subtle coupling be¬ 
tween EVIs and the geometrical constraints involved in 
the comb geometry, this diffusive regime displays a non 
analytical dependence on the vacancy density, checked 
numerically in Fig. 4 of the SM [31]. This is markedly 
different from the case of homogeneous lattices where 
a linear behavior with po is found [32]. In addition, 
the comparison between Eqs. (5) and (11) shows that 
the crossover time between the two regimes behaves like 
tx ~ (Po ln(po))~'*) which can be very large for dense 
systems. As a result, the subdiffusive behavior of the 
first regime is long-lived and potentially observable in 
real systems. We now consider several extensions of these 
results. 

Influence of the initial conditions. The previous results 
were obtained assuming that the vacancies were initially 
uniformly distributed on the lattice. We now assume that 
they are initially located only on the backbone, with a 
linear density Pq" defined as the number of vacancies di¬ 
vided by the length of the backbone. Averaging over this 
initial distribution, which actually amounts to restricting 
the sum over the initial points to the only backbone in 
Eq. (4), it is found that (see SM [31]) 


42n) 


(t) 


lin 

Po 


2V4r(5/4) 




( 12 ) 


Consequently, the time dependence of the cumulants is 
modified in a dramatic way: the cumulants now grow as 
This analytical prediction is successfully confronted 
to numerical simulations (see inset of Fig. 2). This spec¬ 
tacular slowdown of the dynamics with respect to the 
uniform initial conditions is compatible with Eq. (5), 
where now po strictly vanishes. Interestingly, in this case, 
tx —>■ oo, so that there is no cross-over to a diffusive 
regime. 

d-dimensional comb.- The previous results can also 
be generalized to the important case of a d-dimensional 
comb Crf, widely studied in the literature [3, 38, 39], 
which is defined recursively: starting from Ci (a one¬ 
dimensional lattice), Cd is obtained from Cd-i by at¬ 
taching at each point a two-way infinite path (see figure 
in SM [31]). It is found that for uniform initial condi¬ 
tions, the even cumulants all behave like 


lim 

Po-s-O 


K(2")(t) 

Po 


(X t ^ , 
t—^oo 


(13) 


and eventually cross-over to a diffusive linear in time 
regime for d > 2. Note that in the case of the d = 1 
comb, single-file subdiffusion «;(^"^(t) ^ is recovered. 

Finally, reminding that single-file diffusion has been 
shown to be a realization of a fractional Brownian mo¬ 
tion with Hurst exponent 1/4 [40], we conjecture that 
tracer diffusion in a crowded d-comb is more generally 
a realization of a fractional Brownian motion of Hurst 
exponents H = (2^ - l)/2‘^+\ 



FIG. 3. Variance of the TP visiting the teeth, measured for 
different values of po. The collapse of the curves is obtained 
by rescaling time by the crossover time tx,i = l/poi pre¬ 
dicted by our analytical approach, and the variance by p/ 
with 7 ~ 0.62, numerically obtained. The plateau at interme¬ 
diate times confirms the prediction of Eq. (15). For po < 10“^ 
the reported results are obtained via numerical simulations of 
vacancy dynamics, while for po > 10“® via particle dynamics 
(see SM [31]). 


Case of a TP visiting the teeth- We finally come back 
to the original problem of a tracer on a crowded 2-comb, 
where the TP is identical to the bath particles and thus 
allowed to visit the teeth of the comb. The displacement 
of the TP along the backbone can be analyzed in a mean- 
field description that decouples the motion of the TP 
in a tooth from the dynamics of other bath particles as 
a Continuous Time Random Walk, whose waiting time 
distribution ij}{t) describes the time the TP spends on 
a given tooth of the crowded comb. Noting that the 
motion of the TP along a tooth is close to a single-file 
motion, we exp ect that the transverse MSD behaves like 
{y^{t)) cx a/pP [24], in the dense limit. In turn, this leads 
to two different regimes for ip{t): for t tx,i = l/Po> 
{y^{t)) 1, the TP has not had time to explore a tooth 

because of the other crowding bath particles of the tooth, 
and the mean time spent on the tooth is finite; for t 3> 
^x.i, "fit) cx 1/t^, with fi = 7/4, as obtained in [41, 42] 
and checked numerically (see inset of Fig 3). 

The MSD (A^) of the TP along the backbone is then 
related to the MSD of the TP restricted to the back¬ 
bone by using the standard Montroll Weiss relation [34]: 

(^(e) = ^^^3)(^(0). (14) 

Combining the two temporal regimes of K^‘^\t) deter¬ 
mined previously with the two regimes of ^/’(t), we finally 
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obtain that {X‘^{t)) displays three non trivial regimes : 

f if i < ^x,i , 

(X2(t)) OC i t3/4(M-i) =^9/16 ift^_i (15) 

where tx ,2 is a second crossover time (whose explicit de¬ 
pendency on po is not provided by our approach). The 
comparison with numerical simulations shown in Fig. 3 
reveals that: (i) three temporal regimes with expected 
crossover time tx,i are indeed observed; (ii) the expo¬ 
nents of the two first are in good agreement with our 
analytical prediction (15); (iii) the increase of {X^{t)) 


observed in the last regime is in qualitative agreement 
with (15) but the quantitative determination of the cor¬ 
responding exponent would require more extensive sim¬ 
ulations. Remarkably, it is found that in all regimes, the 
dynamics of the TP along the backbone is faster than in 
the absence of bath particles, where {X‘^{t)) ~ In 
other words, the motion of the TP is accelerated along 
the backbone by EVIs. This surprising behavior results 
from two competing effects quantified by our approach: 
hard-core interactions hinder the motion of the TP along 
the backbone but in the same time reduce the time lost 
by the TP in the teeth. 
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Supplemental Material 

I. EVOLUTION RULES OF THE VACANCIES 

In the main text, we first studied the case where the tracer particle (TP) is constrained to move on the backbone 
of the comb. The evolution rules of the particles where given in Fig. 1: (i) the TP jumps on each of the neighboring 
sites of the backbone with probability 1/2; (ii) if a bath particle is located on a tooth of the comb, it jumps on each 
of the neighboring sites with probability 1/2; (hi) if a bath particle is located on the backbone of the structure, it 
jumps on each of the neighboring sites with probability 1/4. The dynamics is constrained with hardcore interactions, 
which means that there is at most one particle per site. 

In the high-density limit, there are very few vacancies on the lattice, and it is more convenient to describe their 
dynamics instead of describing the dynamics of all the particles. The events where two vacancies are on neighboring 
sites are of order O(po^) (where po is the density of vacancies on the lattice). As we focus on the results at leading 
order in po, we do not take these events into account. We adopt a discrete-time evolution, and assume that at each 
time step, each vacancy exchanges its position with one of the neighboring particles. Depending on the position Z of 
the vacancy (see Fig. 4), its jump probabilities are defined as follows: 

• if the vacancy is on a tooth of the comb but not adjacent to the backbone, it exchanges its position with the 
same probability with each of the neighboring particles. 

• if the vacancy is at Z, adjacent to the backbone but not to the TP, the particle located at Z — e-i (resp. Z -|- 62) 
has a probability 1/4 (resp. 1/2) to exchange its position with that of the vacancy. Consequently, the vacancy 
has a probability proportional to 1/4 (resp. proportional to 1/2) to jump in direction —62 (resp. 62), so that 

p(Z-e 2|Z) = Zx i, (16) 

p(Z + e2|Z) = Zx i, (17) 

where p{r\r') is the probability for a vacancy to jump from site r' to site r in a single step, and where Z is 
a normalization constant. With the normalization condition p{Z —)■ Z — 62) + p{Z Z -\- 62) = 1, we find 
Z = 4/3, and finally 


p{Z - e2|Z) = i. 

(18) 

p{Z + e2\Z) = 

(19) 

• if the vacancy is on the backbone but not adjacent to the TP, we can adapt the previous arguments to write 

p(Z±ei|Z) = i, 

(20) 

p(Z±e 2 |Z) = i. 

(21) 

• if the vacancy is on the backbone and at the right of the TP 
vacancy is at the left of the TP), the same arguments lead to 

(this is easily extended to the case where the 

p{Z - ei|Z) = y. 

(22) 

p{Z + ei Z) = -, 

(23) 

p(Z±e2|Z) = 

(24) 

• if the vacancy is on a tooth adjacent to the TP, it steps away 

from the TP with probability 1 as the TP is 


constrained to remain on the backbone. 


These evolution rules are summarized in Fig. 4. 
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II. FIRST REGIME 


In the main text, we established the following expression for the generating function associated with the cumulant 
generating function of the position of the TP: 




po-s-o 


~2po 


— cosk) 

(1-0(1 +A)’ 


(25) 


where we use the simplified notation Fi = F(0|ei;Oi and where +(0|r;0 is the generating function associated with 
Ft{0\r), defined as the probability for a vacancy to reach the origin for the first time at time t knowing that it started 
from site r. We also introduced the quantity i7(0 = Y^^=-oo ^ V 5 ?)■ The determination of the cumulant 

generating function then relies on the determination of H{^) and Fi. Detailed calculations of these quantities are 
given below. 


A. Calculation of H{^) 


The quantity 77(0 is defined as 


= E E ^(0,0 |x,y;O- (26) 

21=1 y — — oo 

The first-passage densities F (0,0|x, 2/;0 are related to the random walk of a vacancy on a lattice with specific and 
non-uniform transition rates, and represented in Fig. 4. 

These first-passage densities can be calculated using the following fundamental relation, resulting from the tree-like 
structure of the two-dimensional comb, and valid for any site r' belonging to the path from r to Tq [43]: 

+(»’ko;0 = F{r\r']OFi'f'\ro]^). (27) 

To compute F (0, Ojx, y; 0 {x > 0), we consider separately the situations where y = 0 and y ^ 0, and obtain: 

• if ?/ ^ 0, the path from {x, y) to (0,0) is decomposed as follows: 


(a;,y) (a:,sgn(y)) (x,0) (1,0) ^ (0,0) 


(28) 










so that generating function of the first-passage density is 

F (0, 0|x, y ; ^) = F (0,0| 1,0 ; ^) F (1, 0\x,0;^)F (x, 0|a:, sgn(y); F (x, sgn(y) |a;, y;^). (29) 

We introduce the quantities 


/ 24 (O = F(l,0|2,0;e), (30) 

/2,2(O = F(l,0|l,l;e), (31) 

/i4(0 = F(l,l|l,2;e). (32) 


The last quantity is relative to the random walk of a vacancy on a tooth far from the backbone, i.e. on a 
one-dimensional structure. It is then given by [34]: 

= (33) 

Finally, using these relations and the fundamental equation (27) once again, we get 

F(0,0|a:,y ;e) = F, [/2.1(e)]'"'-' /2.2(e) [/i.i(e)l'"'-' ■ (34) 


• if y = 0, the path is simply decomposed as: 

(a:,0)^ (1,0)^ (0,0) 


(35) 


and we find 


F (0, Olx, 0 ; e) = F (0,0| 1, 0 ; e) F (1, 0|x, 0 ; e) 

= i^i[/2.i(e)] 


(36) 

(37) 


We now go back to the calculation of 


00 00 


E ^( 0 , 0 |x,y;e) 

x—1 y——oo 


= E 

x—l 


F(0,0|ai,0;e) + 2^F(0,0|x,y;e) 




= FiE/ 2.1(e)"-' 


x—l 

Fi 


i + 2^/2.2(e)/i.i(e)"-' 

/2.2(e) 


1 -I- 2 

i-/2.i(e)V i-/i.i(e)y 


(38) 

(39) 

(40) 

(41) 


We then calculate separately /2.2(e) and /2.1(e): defined in Eqs. ( 30 ) and ( 31 ), which are relative to the random 
walk of a vacancy away from the tracer. 


In what follows, we will use the general relation 

■P(^l’’o;e) = ^^F{r\w;^)p{w\ro), 


(42) 


where p(tu|ro) is the probability to jump from rg to to in a single step, and where the sum runs over all lattice sites. 
This relation is obtained by partitioning over the first step of the walk. 
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1. Calculation of f 2 ,2(0 


We use the general relation (42), and write 

f2,2{0 = ^'^Fiei\w;Op{w\ei + £ 2 ) 

W 

= ^ F(ei|ei;^)p(ei|ei + 62 ) + F(ei|ei + 2 e 2 ;^)p(ei + 2 e 2 |ei + 62 ) 

Because of the tree structure of the lattice, one writes 

P(ei|ei;^) = 1 

P(ei|ei +2e2;0 = ^(ei|ei + e2;0-P(ei + eajei +2e2;0 = /2,2(C)/i.i(C)- 

Finally, we get the following equation for / 2 , 2 (C)' 


/2.2(C) = C 


2 + g/i, 1 ( 0 / 2 , 2(0 


Solving it and recalling the expression of /i,i(C) (Eq. (33)), we obtain 

c 


/2,2(C) = 


I + 2 OI-O' 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


2. Calculation of f 2,1 {^) 


Using again Eq. (42), we write 
/ 2 ,l (0 


(49) 


= C F(ei|ei ;C)p(ei|2ei) + F (ei |3ei; C)p(3ei |2ei) + 2F(ei|2ei + 62 ;C) 7 '( 2 ei + e 2 | 2 ei) (50) 


= C 


g + g/2,l(0^ + 2/2,2(0/2,l (0 


(51) 


/ 2 .i (0 is then the solution of a second-order equation. Choosing the solution fulfilling the condition / 2 ,i( 0 ) = 0, we 
get 


/2.1(C) = ^ - 2 / 2 , 2 (C) - J 2/2,2(C) ) - 1 - 


(52) 


B. Calculation of Fi 


In order to calculate the generating function Fi = E(0|ei;C), we consider the first passage density Ft( 0 |ei) parti¬ 
tioned over the first time step 

Et(0|ei) = ^St,i + = X! Et'_i(ei|2ei)Et_t/(0|ei) + 2 x | ^ Ft'-i{ei\ei -f e 2 )Ft-t'{e\ei). (53) 


t'=i 


i' = l 


The first term of right-hand side corresponds to the situation where the vacancy jumps onto the origin of the lattice 
at the first step (i.e. in direction —ei). The second (resp. third) term corresponds to the situation where the first 
step of the vacancy is in direction Fei (resp. ± 62 ). In terms of generating functions, we get 


El = 


7? 


1 - iCE(ei|2ei ;C) - |CE(ei|ei -f 62 ;C) 

1 -K/ 2 .i(C)-|C/ 2 . 2 (C)' 


(54) 


( 55 ) 
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C. Long-time expansion of the cnmnlants 

Recalling the expression of from Eq. (25): 


we write its expansion in powers of k and obtain 


H{0 


po-o + ^ (2n)!' 


( 56 ) 


(57) 


From the definition of the cumulant generating function '^’{k]^), the coefficients of its expansion in powers of k are 
related to the cumulants of Xt as follows: 


n\ 

n—O 


(58) 


Identifying the expansions (57) and (58), we obtain that all the even cumulants are equal and have the following 
expression 


;(2n) 


iO 


^PoHiO 


m (i_^)(i + Fi) 


(59) 


We also verify that, for symmetry reasons, the odd cumulants are null. Using the explicit expressions of H{^), Eq. (41), 
and Fi, Eq. (55), we obtain the expansion of the even cumulants in the long-time limit (^ —)■ 1): 


-(2n 

hm - ^ 

Po->0 Pq {->1 



1 

25/4(1 _^)7/4- 


(60) 


Using a Tauberian theorem, we retrieve the long-time limit of these cumulants, corresponding to Eq. (5) from the 
main text. 


lim 

PO-i-O 


PO 


~ ^3/4 

t^oo 25/4r(7/4) 


(61) 


D. Influence of the initial conditions 

In the previous sections, we considered the situation where the vacancies are initially uniformly distributed on the 
lattice. The situation where the vacancies are initially located on the backbone of the structure can be studied with 
the formalism presented above. It is straightforward to show that, with this initial condition, the cumulant generating 
function becomes 


where 


^ik;0 


Po->0 


i7'(0(l-cosfc) 

(l_^)(l + i?i) 


iJ'(e) = ^E(0|x,y = 0;0. 

X — 1 

The expression of the even cumulants is deduced as in section IIC, and we obtain 


(62) 


(63) 


j( 2 n) 


(0 


2poFi 


Po-O (l_^)(l + i7^)(l_/2,i(e)] 


(64) 



11 



FIG. 5. Representation of the three-dimensional comb lattice. 


Finally, recalling the expressions of Fi (Eq. (55)) and H'{^) (Eq. (63)) in terms of the first-passage time densities 
computed in section II A, we obtain the long-time limit (^ —>■ 1) of the even cumulants: 


lim 

Po-s-O 


K 


(2n)(^) 


1 


PO 2V4(1_^)5/4- 

Using a Tauberian theorem, we retrieve the long-time limit of these cumulants: 

«;(2")(t) 1 


lim 

Po-s-O 


Po 


2U4r(5/4) 


U/4. 


(65) 


( 66 ) 


E. Definition of a d-dimensional comb 

In the main text, we present the expression of the variance of the position of the TP on a generalized d-dimensional 
comb. This structure is defined recursively as follows: starting from Ci (a one-dimensional lattice), Cd is obtained 
from Cd-i by attaching at each point a two-way infinite path. We represent in Fig. 5 the three-dimensional comb. 


III. SECOND REGIME 
A. Introduction 

The key point underlying the crossover to a second regime is that the previous analytical results are derived by 
taking the limit po 0 before the long time limit. In this case, to leading order, the TP does not move before being 
reached by a given vacancy. In fact, the TP actually diffuses due to its interactions with the other vacancies. Due 
to the effective diffusion of the TP, each vacancy is then assumed to have an additional jump probability in the x 
direction, denoted by D{pq), and which vanishes when po 0. We will write for simplicity e = D{pq). The vacancies 
then perform two-dimensional random walks. We consider simplified evolution rules of the vacancies, which are 
expected to give a correct qualitative description of the system. The evolution rules for the vacancies are represented 
in Fig. 6. 


B. Single-vacancy problem 

As in the previous calculation, we first consider the situation where there is only one vacancy on the lattice. We 
denote by Py\x\Yo-, e) the probability to find the TP at position X at time t knowing that the vacancy started 
from site lEo- This probability can be written in terms of the first-passage time densities Ft(0|r;e) (probability for 
a vacancy to reach the origin for the first time at time t starting from r) and the conditional first-passage densities 
F)*(0|e^|r; e) (probability for a vacancy to reach the origin for the first time at time t starting from r and being at 
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FIG. 6. Simplified evolution rules of the vacancies, when they experience an additional jump probability e = D{po). 


site at time t — 1). Summing over the number of steps p taken by the TP, over the directions z/i,..., i/p of these p 
steps and over the length of the time intervals mj elapsed between the (j — l)-th and j-th step, we get 


Pi^\x\Yo;e)=Sx,o 1 - ^F,(0|ro;e) 


i=o 

OO OO OO OO 

+ EE-E E dmi+...+mp+i,t E--E ^e„j^+...+epp,X 

p=l ?7ii = l mp — lmp+i—0 1^1 I'p 

( rup+i \ 

1 - 


The Fourier-Laplace transform of Eq. (67) writes 


^(1) 1 
P ik\Z;^,e) = 


1-^ 


1 + 




(67) 


( 68 ) 


where we defined: 


= det[l - r(fc;^,e)] (69) 

( e''=P*(0|ei|e_i;e,e) e''=P*(0|ei|ei;C, e) \ 

T(fc;5,e)= (70) 

^e-'^^P* (0|e_i|e_i; e) e-i'=P*(0|e_i|ei; C, e )) 

U,{k; e, e) = V{k- e, e) E (1 - a,6 (^1) 

a—ibl 

The TP being symmetric, we have the following relations 

P*(0|ei|ei;e,e)=P*(0|e_i|e_i;e,e) (72) 

P*(0|ei|e_i;e,e)=P*(0|e_i|ei;e,e), (73) 

so that 'D{k;^,e) and P±i(fc;^,e) reduce to: 

P(fc; e, e) = 1 + P*(0|ei|e_i; e, e)^ - P*(0|ei|ei; ef - 2P*(0|ei|e_i; C, e) cosfc, (74) 

U±i{k- e, e) = (e^**^ - 1)[1 - eT*'=P*(0|ei|e_i; e) - F* (0|ei |ei; e)]. (75) 
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We finally obtain the following expression of the single vacancy propagator in terms of the first-passage time densities 
F*: 

P ik\Z-^,e) = — 1 +- - ,2 .2 ^ ^ - - > (76) 

l-f-F'n +[Fli) -2Fl^cosk 


where we introduced the simplified notation 


F; =F*(0|ei|e,;^,e). 


C. Finite density of vacancies 


As for the study of the hrst regime, we now assume that the density of vacancies is equal to po, and we show that 
the Fourier-Laplace transform of the cumulant generating function in the limit po 0 is given by 

^ Po ^ ^F*(0|eJZ;e,e). (78) 


The TP being symmetric, the quantity ?’^) independent of a, 

^ F*(0|e_i|Z;^,e) = ^ F*(0|ei|Z; C, e), 


and we obtain 


~^-po (5]F*(0|ei|Z;C,e) _ Pe)e- 

\Z#0 / a=±l L ? 


Defining S(^, e) = X)z 3 io finally have 


$(fc;e,e) ~ -poS(e,e) ^ ^ _ P'\fc|e_a;e)e- 

Po ->0 [1 ~ S 


The expression of the variance of the position of the TP is easily deduced from 


K^^TC>e) = - 


a2$(fc;^,e) 


Using Eq. (76), we get 


«(2)(^,e) = _2poS(e,e)- 


F* _ p*^ _ 1 


- 1 -f Fli [F* + 1- Fl^ 


Consequently, the determination of the variance of the TP position only relies on the estimation on the following 
quantities: 

• the sum S(^,e) = Ez#o(0 |ei |^; C. e), 

• the first-passage time densities F^i = F*(0|ei|e±i; e). 


In what follows, we evaluate the Hrst-passage time densities F^i as well as the sum E(^, e). 
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1. Evaluation of F*{0\ei\e.±\-,^,e) 

In what follows, we write the conditional first-passage densities F* (0|ei |e±i; e) in terms of the propagators of the 
random walk of a vacancy with the rules specified in Fig. 6. We first aim at computing (0|ei|e_i; e). We then 

assume that the vacancy starts at site e_i, and we write the following partition over the visits of the vacancy to the 
origin of the lattice: 


t-i 


p{ei —)■ 0)Pi_i(ei|e_i; e) — P/(0|ei|e_i; e) -I- p(ei —> 0)Pi/(0|e_i; e)P(_i_t/(e_i|0; e) 


(84) 


t'=0 


where p(r —>■ r') is the probability for a vacancy to jump from site r to site r' in a single step (here p{ei —>■ 0) = l/4-|-e). 
Multiplying the previous equation by and summing for t going from 0 to oo, we obtain the following relation between 
the associated generating functions: 


C +e) ^(ei|e_i;^,e) = P*(0|ei|e_i; e) -h ^ 

The first-passage density F(0|e_i;^, e) is evaluated with the following renewal equation [34]: 

P(0|0;f,e) 

We obtain the following expression of P*(0|ei|e_i; e) in terms of the propagators P: 


P*(0|ei|e_i;^,e)=^( -+e 


P(ei|e_i;^,e) - 


P(0|e-i;^,e)P(ei|0;g,£) 

P(0|0;$,e) 


(85) 


( 86 ) 


(87) 


The random walk performed by the vacancy being translationally invariant in the direction of the backbone, we will 
use the following relations: 

P(ei|e_i;e,e)=P(2ei|0;e,e) (88) 

P(0|e_i;e,e)=P(ei|0;^,e). (89) 

We define a{r;^,e) as 

air; e) = P(0|0; e) - P(r|0; e). (90) 

Then, the conditional first-passage densities are simply expressed in terms of Q;(ei;^, e), a(2ei;^, e) and P(0|0;^,e): 


P*(0|ei|e_i;^, e) — ^ ( - -be 


2a(ei;^,e) - Q;(ei;^,e) - 


a(ei;^,e)2 


P(0|0;e,e) 

A similar calculation leads to the following expression of P*(0|ei jei; 5, e) in terms of a(ei;^,e) and P(0|0;5,e): 


(91) 


P*(0|ei|ei; e) — ^ ( - -be 


2 a(ei;^,e) - 


«(ei;g,e)^ 

P(0|0;e,e) 


(92) 


In what follows, we compute the quantities a(ei; e), Q;(2ei; e) and P(0|0; e) associated with the random walk 
of a vacancy with the evolution rules presented in Fig. 6. We follow the method first introduced by Nieuwenhuizen 
et al. [36]. 

For generality, we assume that if the vacancy is on the backbone (respectively not on the backbone), it has a 
probability p{j (resp. pH) to go left or right, and (resp. p^) to go up or down. Let Px,y{t) be the probability to 
find the vacancy at site (x, y) at time t, knowing that it started from the origin. The master equations of the problem 
are the following: 

• if 2/ 7^ 0,±1: 


Px,yit + 1 ) — P^^ [Px+l,y{t) + Px-l,yii)\ + P'^[Px,y+l{t) + Px,y-l{i)\ 


( 93 ) 
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• if y = 0: 


Px,o{t + 1) — P'^[Px,i{t) + Px,-iit)] + pf [Px+i,oit) + Px-i,o{t)] 


(94) 


• a y = ±1: 


Px,±i{t + 1) 


P'^ Px,±2it) + p-t Px,oit) + p^^[Px+i,±iit) + Px-l,±l{t)]. 


(95) 


We introduce the following generating functions and Laplace transforms: 

OO OO 

p,{kx-,o = T. E 

t —0 x— — (X) 

OO OO 

P{kx,ky-,0 = J2 E 

t=0 x^y——oo 

Using the master equations, it is straightforward to show that P and Pt are related by 


P{kx,ky]C) 


1 + -p^)cosky + {pI -p\'') cos kx]Pb{kx-,0 
1 — 2^[p-*- cos ky + pll cos kx] 


(96) 

(97) 


(98) 


In order to get an equation for Pi,(kx',^), we integrate each side of this equation over ky and use the simple relation 
between P and Pf,- 


/•27r jjL 

Pb{kx;0 = J ^P{kx,ky-,0- 


We also use the following useful integrals: 


dk 


/q 27r cosh /i — cos k sinh p ’ 


dk cos k 




/q 27r cosh p — cos k sinh p ’ 


to finally obtain 


Pb{kx;0 = 


2^p-^ sinh/X — + {pi — ;xll) cos 


where p is such that 


, 1 — 2^pll cos kx 

= -2iP-^ 

In what follows, we will use the explicit expressions of the transition probabilities: 

pll = e, 




Pt = 4 - e- 

a. Computation o/a(2ei;^,e) We focus on the computation of Q;(2ei;^,e), defined by: 

a(2ei; e, e) = P(0|0; C, e) - P(2ei|0; e, e) 

('2tt 




(99) 

( 100 ) 

( 101 ) 

( 102 ) 


(103) 

(104) 

(105) 

(106) 

(107) 

(108) 
(109) 
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In the limit where we first take ^ = 1, one has: 

PbiK; e = 1) =-=- ~ ^ 

1 — coskx + 2-\/e(l — 4e)-\/(l — cos A:a;)(l — e — ecos kx) 

and consequently 

^ 2(l-2e)(l-e-2i'=^) 

a(zei; 4 = 1, e) = / —-z- , _. 

Jo 27r 1 — cos fca; + 2-\/e(l — 4e)-\/(l — cos fca;)(l — e — ecos fca;) 

For symmetry reasons, the imaginary part of the integral is zero, and 

.O . 1 ^ ^ r (l-2e)sin^fc, 

a(2ei;4 = l,e) = — / akx - , - . 

Jo 1 — cos kx + 2-v/e(l — 4e)-\/(l — cos A:a;)(l — e — ecoskx) 

With the change of variable u = cos kx, one gets 


a(2ei;C = l,e) = 


a( 2 ei;^ = l,e) = - / dkx 
7r ,/n 


4 Vl — 

.(2e.;| = 1,0 = -d - 20/__ '^" 1 _ „ ^ _ „)(i _ , _ 

= -0-2e) f du - ^ “ 

TT J_i -y/l — M + 2-v/e(l — 4e)\/l — e — eu 


= ^(l-2e)/ dw^^- 7^"" ^_ _ . (114) 

TT J_i -y/l — M + 2-y/e(l ~ 4e)\/l — e — eu 

We aim to compute the expansion of a(2ei|0;^ = l,e) in powers of e. The integral in the previous equation have 
an explicit expression: 

a(2ei;^ = l,e) =- ~ -^[256v^7(l - 2e)e9 - 64^27(1 - 2e)e7l6v^7(l - 2e)e5 + 327Te^ 

7r(2e - 1) (32e3 + 2e + 1)" 

+ 4-\/2a/( 1 — 2e)e3 + 2i-\/27r7l6e^ + e — 327re^ — 2(1 — 4e)^ (8e^ + 2e + l) arctan 
— 8i-\/27r7l6e^ + e^ + 4-\/2 f 7l6e^ + e — 4\/l6eM-^) arctanh [ ) 

V /V d/-.'12e3 + 16e2-2e + l / 


/-32e3 + 16e2 - 2e + 1 


+ 27re-2v^7(l-2e)e-7r]. 


Its expansion for e —)■ 0 leads to: 

o /o 

a(2ei;^ = l,e) = 4H- •\/elne+--- 

£->0 TT 

b. Computation o/a(ei|0;^ = l,e) Similarly, we study a(ei|0;,J = l,e), defined by 

a(ei;^,e) = P(0|0 ;C,e) - P(ei|0 ;^,e) 

Using again Eq. (110), and with the change of variable u = coskx, we obtain 


a(ei;Cu) =-(1 - 2e) / du— 

TT J_i -y/l 


- m2 + 27(1 - 4e)7(l + m)(1 - e - eu) 


The explicit expression of this integral is 


= 01 + 40(1 - 2. + 8i0Vi'/l + ^ -- + !) + 4i«(4. - 1) 

+ Arcsin(4e — l)7l + 16e22-\/2e^/^(4e — 1) + 2e(l — 4e) Ine + 4e(l — 4e) ln(l — 4e) 
+ 2e(4e - 1) ln(-16e3 + + 1 + 86^ - e + 1)]. 

Its expansion for e —>■ 0 leads to: 


Q;(ei;^ = l,e) = 2-\ - ^/l\ae + - 

£-40 TT 


( 121 ) 
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c. Computation o/P(0|0; e) Using the definition of Pb from Eq. (96), the propagator P(0|0; e) can be written 

dK 


P{0\0;^,e)= -^Pbik^;a 

where Pb{kx',0 is defined by Eqs. (102) to (107). We then express P(0|0;^,e) as the following integral 

dkx _ 2(1 - 26 ) _ 

/q 27r 1 — ^coskx + (1 — 4e)-\/(l — 2^ecoskxY ~ ~ ^e)^ 

With the change of variable u = cos kx, this integral is rewritten as 


P(0|0;^,e) = 


P(0|0;^,e) = - 


du 


1 - 2e 


ttVo Vl - 1 - + (1 - 4e) v^(l - 2^eu)‘^ - ^^(1 - 2e)2 ’ 


( 122 ) 


(123) 


(124) 


In the limit where ^ 1 and e —)■ 0 this quantity diverges. Studying the integrand, we show that this divergence 

is located near u = 1. Introducing a new integration variable v = 1 — u, we obtain that, for ^ 1, the integral is 

estimated at leading order by 

r-2 


P(0|0;e,e) 


1 - 2e 


dv 


5^1 7re(l 4e) Jq ^^2 - v^v + ^ “ 2 

The integral over v is conveniently expressed in terms of the complete elliptic integral of the first kind K(cc): 


P(0|0;^,e) 


5^1 7re(l - 4e) ^1 - ^ 
We use the following useful expansion of K(x): 


l-2e 2^/2e ^(2e{Ae-l-C) 


1-? 


K - 


1 


, ~ ^Inl, 

X) a;-i.0+ 2 X 

to finally obtain the expansion of P(0|0; e) at leading order in (1 — ^): 

?(0|0;f,e) 




4-i-i •K\f2t(\ — 4e) 1 ~ C 


(125) 


(126) 


(127) 


(128) 


Finally, recalling the expressions of P*(0|ei|e±i; e) in terms of the functions a and of the propagator P(0|0; e) 
(Eqs. (91) and (92)), we use the results from Eqs. (121), (116) and (128), we obtain the following expansions 


1 


P*(0|ei|e±i; 5, e) — ^±i(e) + C{e)- -j-h ..., 

?->i In 


1 -? 


where 7lj-i(e) and ^(e) have the following expansions when e —>■ 0: 

^i(e) 1 


e-J-O 


A_i(e) ~ — v^ln- 

e-s-O TT e 

C{e) ~ -TTv^. 

€—>•0 


(129) 

(130) 

(131) 

(132) 


Recalling that the additional jump probability e = P(po) is a vanishing function of p when rhoo —>■ 0, the conditional 
first-passage densities can be rewritten in terms of pq as follows: 


P*(0|ei|e±i; po) — ^±i(po) + Ct{po)- -j-h ..., 

In ^ 

where kl±i(po) and C'(po) have the following expansions when po —>■ 0: 

^i(Po) ~ 1 

€—>■0 

A_i{po) oc — y/D{po)ln—^ 

e^O TT D(Poj 

C{po) oc -^^^yD{po). 

€—>■0 


(133) 

(134) 

(135) 

(136) 
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H. Evaluation of 


The evaluation of the sum S(^,e) is conveniently done by adopting a continuous-space description of the system. 
Following the approach proposed by Arkhincheev and Baskin [5, 44], we dehne by G{x,y\xo,yQ]t) the propagator 
associated with the random walk of a vacancy starting from site (xq, j/o) and arriving at {x, y). The diffusion equation 
verified by G is then 


— - [D,5{y) + D{p,)] — , 

^G{x,y\xo,yo;t = 0) = 6{x - xo)6{y - yo), 


(137) 


where Di (resp. D 2 ) is the diffusion coefficient of the vacancy in the x (resp. y) direction, and D{pq) is the additional 
jump probability experienced by the vacancy due to the effective diffusion of the TP. This coefficient is not explicitly 
known, and is only assumed to vanish in the limit po —>■ 0 . 

We now solve the diffusion equation (137) in order to obtain an expression for G{x,y\xo,yo;t). We introduce the 
continuous Laplace transform, defined for any time-dependent function by 

7(s) = / dte~^*f{t), (138) 

and the Fourier transform along the x-direction, defined for any x-dependent function: 

/ oo 

dxe"'''^ 5 (x). (139) 

-00 

The Fourier-Laplace transform of (137) is 

G{k,y\xo,yo;s) = e~^'^^°6{y - yo). (140) 


s + D{po)k^ — 


We first assume that yo > 0. On each of the three domains ] — 00 , 0[, ]0,7/o[ and ]j/o, cx)[, Eq- (140) becomes 

52 T 


s + D{po)k -£> 2 ^ 


G{k,y\xo,yo;s) = 0, 


(141) 


whose general solution is 


I \ .1 I s + D{po)k‘^\ I ls + D{po)k'^ 

Gik,y\xo,yo;s) = Aexp \-J -- -ySexp' 


D2 


D2 


(142) 


where A and B are two constants to be determined. The physically relevant solutions on the three domains take the 
form 


G{k,y\xo,yo\s) 


< ^2 exp + B 2 exp 


if y < 0 , 
if 0 < y < 2 / 0 , 
if 2/0 < y. 


(143) 


There are four constants (i?i, A 2 , B 2 and A 3 ) to be determined. The continuity of G as a function of y yields the 
following relation: 


Bi = A2 + B2, 

Ao = A 2 + B 2 exp , 


( 144 ) 
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so that we write 


G{k,y\xo,yo-, s) = <( 


A 2 exp + B 2 exp 


l s+D(po)k'^ 

D 2 y 


'Ay < 0 , 

ifO<y<yo, (145) 


A2 exp + -^2 exp \ - y)^ if yo < y. 


The constants A2 and B2 are obtained by integrating Eq. (141) respectively over the intervals [—rj, rj] and [yo — V^ Vo + v] 
and by taking the limits 77 —)■ 0. On the intervals ] — 00 , 0[ and ]0, yo[, the following relation holds 


G(fc,j/|xo,j/o;s) = 2l2exp 


I s + D{pQ)k‘^ 
142 


\y\ +B 2 exp 


I s + D{po)k‘^ 
^2 ' 


We integrate Eq. (141) over [— 77 , 77 ] and using the following general relation (which holds for any a > 0) 


r A r 


^-Ay\ 


dy a[a — 25{'y)]e —>■ —2a 


7]—^0 


we obtain the relation 


Dik^(A2 + B 2 ) + 2A2^/D2[s + D{po)k‘^] = 0. 

On the intervals ]0,7/o[ and [j/o,cx)[, the following expression of G holds: 


(146) 


(147) 


(148) 


G(fc,2/|a;o,J/o; s) = ^2exp 


I s + D{po)k‘^ 

^2 


y +B 2 exp -A 


I s + D{po)k‘^ 


D 2 


-yo exp 


I s + D{pQ)kP- 

^2 


We integrate Eq. (141) over [7/0 — '' 7 , 2/0 + v], and using again Eq. (147), we obtain the expression of B2: 

1 


\y-yo\j ■ 

(149) 


B 2 = 


2^D2[s + D{po)k^] 

The relation between A 2 and B 2 yields 

Dik^ 


exp —ikxo — 


' s + D{po)k‘^ 


-yo 


(150) 


Ao = — 


Dik^ + 2y/D2[s + D{po)k^] 2y/D2[s + D{po)k^] 


exp —ikxo — 


I s + D{po)k'^ 


Do 


-yo 


(151) 


In what follows, we will only consider the propagator of a random walk starting from a generic point (xo, 2 /o) and 
arriving at the origin (0,0), so that we will take y = 0. Generalizing the calculation to the case of a starting point 
located in the domain j/g < 0, we finally obtain the general expression 


G(fc, 0 |xo, 2 /o;s) = 


exp (^-ikxo - 

Dik"^ + 2y/D2[s + D{po)k‘^] 


(152) 


We now study the asymptotic behavior of the sum E(^,e) = in the limit where ^ 1 

(long-time limit). In the continuous-space and time description, we replace the sum by the following integral 


S(s,Po)= / dxg / d 2 /oF( 0 , 0 |a;o, 2 /o;s)- 


Writing the following renewal equation 


E(0,0|a:o,2/o;s) = 


G{0,0\xo,yo;s) 

G(0,0|0,0;s) 


(153) 


(154) 
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we get 


S(s,Po) = 


G(0,0|0,0;s) 

Taking the inverse Fourier transform of Eq. (152), we write 


dxo / di/oG(0,0|xo,j/o;s)- 


(155) 


^ /•o° ^1^ 

G{x,0\xo,yo;s) = —exp{ikx) 


exp ^-ifcxo - 
£»ifc2 + 2y/D2[s + D{po)k^] 


I 2 / 0 I 


(156) 


We now evaluate separately G(0, 0|0,0; s) and / dxo dyo G(0,0|a:o, Vo] s): 
• using Eq. (156), we obtain the expression of G(0,0|0, 0; s): 


G(0,0|0,0;s) = 


dk 1 

^ D^k^ + 2y/D2[s + Dipo)k^] ’ 

- f dk - ^ =. 

TT Jo Dik^ + 2^D 2 [s + D{po)k^] 


(157) 

(158) 


In the long-time limit (s —>■ 0), the integral is dominated by the small values of the variable k. Introducing a 
cutoff value A, we write 


G(0,0|0,0; s) ~ — / dfc 


1 


ttJo 2^D2[s + Dipo)k^]' 


Using the following integral (for a > 0), 

(\x 


\/a? +: 


= In 


A A 


a \ a 


+1 , 


we obtain the leading order behavior of G(0, 0|0,0; s) when s —)■ 0: 

1 


G(0,0|0,0;s) 


^ In-. 


s-i-o dn 1 / 1)2 \JD[po) s 
the integral f dxo / dyo G(0,0|a:o, 2 /o; s), using Eq. (156), yields: 

^ exp (^-ikxo - 


dxo / dyoG( 0 , 0 |a:o, 2 /o;s) = / dxo / dyo 


dk 


|j/o| 


r°° dk /• , 
— / dxo 


27r Bik^ + 2^D2ls + D(po)k2] 
exp(—ifcxo) I 


Do 


-00 ^ J Dik^ + 2^D2 [s + Dipo)k^] V s + D{po)k^ 
2Sik) I 


dk 


Do 


Dik^ + 2yJD2 [s + D{po)k^]\ s + D{po)k^ 


Finally, using the estimates from Eqs. (161) and (165) in Eq. (155), we get 

Atti/DJi/d/po) 

Ss,po ~ -wTTT-■ 

s-j-o sln(l/s) 

In the long-time limit, we get the equivalent discrete time description in terms of the variable /. 

d'n\jD 2 \J D{po) 




(l-Olny^ 


(159) 


(160) 


(161) 


(162) 

(163) 

(164) 

(165) 

(166) 

( 167 ) 
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FIG. 7. (Color online) Variance of the TP when its motion is constrained on the backbone, rescaled by the time t and 
Po log(l/po)^, according to Eq. (11) of the main text, for several values of po- Inset: the black dots represent the diffusion 
coefficient D{po) = /2t measured in simulations and the black line is the best fit obtained with the function 

f{x) = ax^log{l/x)^ (a=0.2336). For comparison we also show the best £t obtained with the function g{x) = bx^, with 
b = 1.375 (red line). 


3. Expression of the variance 

We recall the expression of the generating function associated with the variance of the TP position: 

= -2poS(C,po)- 


Ff - - 1 


[F* - 1 + Ffi 




1-F 


-1 


(168) 


Using the results from Eqs. (167) and (133), we obtain the expansion of the generating function associated with the 
variance up to a numerical prefactor, where we first take the long-time limit (^ —i 1) and ultimately the high-density 
limit {po —)• 0); 


1 




(169) 


Using a Tauberian theorem, we retrieve the time-dependence of the cumulant and the expression presented in the 
main text: 


lim (X po\/D{po) In 

t^oo t Po—>-0 


(170) 


From this expression we obtain Eq. (10) of the main text, which turns to be in good agreement with numerical 
simulations, as shown in Fig. 7. 


IV. NUMERICAL SIMULATIONS 

For the case where the TP is constrained on the backbone, the data reported are obtained via Monte Carlo numerical 
simulations of vacancy dynamics. We considered a lattice of size x Ly = 1000^ with periodic boundary conditions 
and a number of vacancies M = poL^Ly. At each Monte Carlo step, the position of all vacancies is updated according 
to the evolution rules described in Section I. If a vacancy exchanges its position with the TP, the TP position is also 
updated and the displacement is measured. Reported data are averaged over some thousands of realizations. 

In the case where the TP is allowed to visit the teeth of the lattice, for values of po < 10“^ the reported data are 
obtained via Monte Carlo simulations of vacancy dynamics, with the same lattice size as above. For larger values 
of Po we performed Monte Carlo simulations of particle dynamics. At each Monte Carlo step, all particles (included 
the TP) attempt a move on a neighbor site with the probabilities given in the main text, and the move is accepted 
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if the target site is empty. In this case, the lattice size is x Ly = 200^ with periodic boundary conditions and the 
number of particles is fV = (1 — po)LxLy. 



